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1. Introduction 

In this paper we propose a construction of a first order parametrix for solutions to 
the covariant, tensorial wave equation 

□g* = F, (1) 

under minimum assumptions for the Lorentz manifold (M, g). Here ^' and F are 
k tensor- fields on a 3 + 1 dimensional Lorentz manifold (M,g) and 

denotes the covariant wave operator on M, with D the Levi-Cevita connection 
defined by g. To simplify the discussion below we consider first the scalar case 

□gV- = /■ (2) 

In Minkowski space (M^+^, m) with m = diag{ — 1, 1, . . . .1} the wave operator on 
the left hand side of (2) is the standard D'Alembertian □ = m°'^dadi3. The general 
solution of Oip = f can be written in the form, 

V' = V/ + V'o (3) 

with ipa a solution of the homogeneous equation Dipo = and ipf given by the the 
Kirchoff formula, 

i;fit,x) = (47r)-i / / \x-yr'fis,y)dsdaiy) 

Jo J\x~y\=t-s 

= (47r)-i / ^—S{t-s-\x-y\)f{s,y)dsdy. (4) 

Jm.\+^ f - y\ 

Here da{y) denotes the area element of the sphere — ?/| = t — s and S represents the 
one dimensional Dirac measure supported at the origin. The homogeneous solution 
ipo is fixed by initial data on the hyperplane t — 0. 

One can also recast (4) in the form 

V'/(t, x) = (Att)-' f Hit - 3)d{ -{t- + \x- y)dsdy (5) 
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where H{t) is the Heavyside function supported on the positive real axis and the 
expression |a; — — [t " s)"^ = <io{p,q)'^ is the square of the Minkowski distance 
function between the vertex p = {t,x) and the point q — {s,y) in the causal past 
J~{p) n Ri^'^^ of the point p G R"^^^. All attempts to extend Kirchoff's formula to 
a general four dimensional curved space-time are based on either (4) or (5). Thus 
the first term in the so called Hadamard parametrix is constructed by replacing 
the Minkowski distance function do with the Lorentzian distance function d{p, q) 
defined by the metric g. Thus one can set, 

V-zb) = (4^)-^ / rip, q) 5{dHp, g)) /(g) d;v{ji) (6) 
-iJ-ip) 

with dv the volume element of the metric g, and r{p, q) a correction factor which 
verifies a transport equation along the null boundary of J^{p) and such that 
riPjP) = 1- The integral on the right makes sense for the portion of J'^{p) which 
belongs to a neighborhood V oip where the geodesic distance function d{p, q) is well 
defined and sufficiently smooth. Typically one requires V to be causally geodesi- 
cally convex, i.e. any two causally separated points in V can be joined by a unique 
geodesic in V. The local parametrix in V is then defined 

^/(p) = (47r)-i / r(p, q) 6{d\p, q)) f{q) dv{q) (7) 

Jj-ip)nv 

The integral in (7) is supported on the portion of the boundary J\f~ (p) of J'~ (p) 
included in T>. 

The error term O^tpf — f, however, does not vanish unless g is the flat metric m. 
One can improve (6) by making successive corrections based on solving a series 
of transport equations in J'^ (p) n T>. In the process the error term can be made 
as smooth as we wish, for given regularity of /, at the price of requiring higher 
regularity of the metric g, see [Pried]. Moreover the resulting parametrix, called 
Hadamard parametrix, is no longer supported just on the boundary of J'~{p). One 
obtains a solution of (2) of the form, 

V'(P)= / E_ip,q)f{q)dv{q). (8) 

Jj-{p)n-D 

with E-{p, q) = r{p, q)5(^d?{p, g)) + . . . is the retarded Green function of Dg 

The Hadamard parametrix (8), which requires both infinite smoothness of g and 
geodesic convexity for V is ill suited for applications to nonlinear problems. It 
turns out that in many situations one does not need the precise representation 
(8) and that in fact the first order parametrix of type (6) sufiices. This fact was 
first made use of by S. Sobolev, see [Sob], to provide a proof of wcU-poscdness 
for general second order linear wave equations with variable coefficients. A similar 
parametrix was later used by Y. C. Bruhat, see [Br], in her famous local existence 
result for the Einstein vacuum equations. Both [Sob] and [Br] construct their first 
order parametrices, which we refer to as Kirchoff-Sobolev, based on the flat space 
formula^ (4). The generalization of (4) to a curved space-time proceeds from the 



It is easy to show that the two constructions (6) and (13) differ in fact only by a normalization 
iaictor at the vertex p. 
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observation that the function Up(s, y) = t — s — \x — y\ is an optical function, i.e. 

ni°'^daudf3u = 0, (9) 

vanishing precisely on the past null cone J\f~ (j>) with vertex at p = {t, x) given by 
the equation Up = 0. Letting r = |a; — j/| one can easily check that 

D{r-^5{up)) = {Dr-^)S{up) + {-2L{r-^)+r-^Dup)6'{u) 

+ {ui''^docUpd0Up)5"{u) = 47r5{p), 

with S{p) the four dimensional Dirac measure supported at p. Indeed the terms 
involving 5"(up) and S'{up) both vanish, the first in view of (9) and the second 
because, 

-2L(r-^) + r"^ Dup = 0, 

with L the null vcctorfield along M~{p) defined by L = —m^^d^Upda- On the 
other hand d{up)Dr^^ = 6{up)Ar~^ = 47r(5(p). 

Based on this one can generalize (4) to a curved space-time by setting, 

i>f{p)= ( a{p,q)d{up{q))f{q)dv{q) (10) 

JJ-{p)n-D 

where Up = Up{q) is the backward solution to the eikonal equation, 

g"'^a„u9^u = o, (11) 

vanishing on the past null cone 7V^(p), and a{q) = a{p,q) verifies the transport 
equation similar to that satisfied by r^^ in flat space. As in (6) we need to restrict 
ourselves to a neighborhood 2? of p in which solutions to (11) remain smooth. 

To explain the restriction to the neighborhood V to which the integral in (10) is 
restricted we return for a moment to the initial value problem in flat space-time. In 
the Minkowski space-time model with the choice of an initial Cauchy hypersurface 
So = = 0} the Kirchoff formula 

i,fip) = (47r)-i / _±_Siupiq))f{q) dv{q) (12) 

with p ^ {t,x),q = (s, y), Up{q) = t — s — \x — y\ and r{p, q) = \x — y\, coincides at 
point p with the solution of Oip = f with zero initial data at t = 0. The represen- 
tation is valid for any point p to the future of Sq and the surface of integration 

Ar-(p)n J+(So) = {(s,y) : t-s = \x-y\, s>0} 

is smooth with exception of the vertex point p. In a flat space-time model with 
the Lorentzian manifold M = R x 11^, where Ha = x R/aZ is a flat cylin- 
der of "width" a, the representation (12) also coincides with the solution of the 
inhomogencous wave equation with zero initial data at t = 0, provided that we 
restrict ourselves to points p = {t, x) such that t < a. For points p = {t, x) with 
t > a formula (10) no longer^ represents the solution of the inhomogencous problem 
with zero initial data at t = 0. The null hypersurface M~{p) fl J^{Y,o) develops 

^In fact the correct representation can be obtained by lifting the problem to the covering space 
R X R'^ X R, applying the Kirchoff formula and taking periodization in the last variable with the 
period a. 
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singularities'^ (scars) in the time interval [0,i — a] duo to intersecting null geodesies. 
This shows that the accuracy of the Kirchoff formula in this case is restricted to 
the neighborhood V = {{t,x) : < t < a} oi the Cauchy hypersurface Eq. 

To describe the situation in a general sp;u;e-time (M, g) we assume that M is glob- 
ally hyperbolic, i.e., there exists a Cauchy hypersurface S C M with the property 
that each in-extendible past (future) directed causal curve from a point p to the 
future (past) of E intersects E once. We denote by E+ = J7'"'"(E) the future set of 
E. By finite speed of propagation the solution tp{p) of the wave equation d^ip = f 
at point p G E+ is completely determined by the values of / in J~{p) Ci E_|_ and 
initial data for ^ on J~{p) fl E. 

Definition 1.1. We will say that E-{p,q) is the retarded parametrix for Dg at p 
if 

Mp)= I E_{p,q)f{q)dv{q) 

coincides with the solution of the problem Dg^ = / with zero initial data on 
E. We will say that the first term in the expansion of E-{p,q) - distribution 
K~ = a{p, q)5 {up{q)) - is the retarded Kirchoff-Sobolev parametrix. 

Let I? be a space-time neighborhood of E. The expression 

i>}{p)= I a{p,q)5{up{q))f{q)dv{q), p&V (13) 

Jj-{p)nY.+ 

is the Kirchoff-Sobolev approximation to the solution ^(p) of the wave equation 

'^gV' = f with zero initial data on E. Clearly tpj fails to be a solution to (2) in the 
non-flat case. We write a general solution of (2) with zero initial data on E in the 
form, 

^{p) = ^f{p)+£f{p) (14) 

with £f an error term. 
In this paper we will: 

(1) Provide a careful derivation of (13) and (14) for points p in a suitable 
neighborhood P of E and show that the error term £f can be expressed in 
the form, 

£f{p)= ( £{p,q)5{up{q))i^{q)dv{q) (15) 

where the smooth density £{p,q) depends only on geometric quantities 
associated to the null hypersurface J\f~{p). 

We should note that classical constructions of the Kirchoff-Sobolev parametrix 
establish the error term £f as explicitly dependent on the metric g and its 

derivatives relative to some chosen system of coordinates. To our knowledge 
the fact that £f is supported only on the boundary of the past set J~{p) 



Note that although past null geodesies intersecting, say at t, = t — a can be extended beyond 
t* they no longer belong to the boundary of the causal past of p. 
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does not seem to have been fully recognized and used in applications. A 
similar observation was, prior to this work, communicated to \is verbally by 
V. Moncrief. His claim, based on Friedlander's treatment of the Hadamard 
parametrix, was the starting point of our own investigations. 

(2) Extend formulas (13) and (15) to the covariant tensorial wave equation (1). 

Once again the classical treatment of the tensorial wave equation introduces 
additional coordinate dependent error terms. Our approach is entirely co- 
variant. 

(3) Provide a minimum set of conditions for the local geometry of M near p to 
ensure that the representation (13) and (14) holds true at p. We also make 
use of our recent results from [Kl-Ro4] to show that for the Einstein vacuum 
space-times (M, g), with vanishing Ricci curvature, formulas (13) and (14) 
can be extended to points p at distance from S, with dependent, 
essentially, only on the norm of curvature* of g. 

(4) Our formula can be easily adapted to gauge invariant wave equations. In 
section 4 of the paper we write down such a formula and show how it can 
be used to give a very simple proof of the Eardley- Moncrief global existence 
result for the Yang-Mills equation in the 3-1-1 dimensional Minkowski space, 
sex; [EM1],[EM2]. The remarkable fact about our approach is that it is 
entirely gauge independent; we don't need to specify any gauge condition^. 

The size of the neighborhood V, mentioned above, is first and foremost constrained 
by the condition that the optical function u is smooth. In the case of a Riemannian 
manifold the distance function from a point p is smooth in a geodesically convex 
neighborhood of p whose size can be evaluated in terms of the norm of the metric 
g, as measured in a given system of coordinates. Alternatively, by a theorem of 
Cheeger, the size of this neighborhood depends only on the pointwise bounds for 
the Riemann curvature tensor and a lower bound on the volume of a unit geodesic 
ball. For similar reasons the construction of a solution Up to (11) is restricted to a 
geodesically convex neighborhood^ T> of p. Unlike the Riemannian case, however, a 
purely geometric characterization of the size of a geodesically convex neighborhood 
of a point p is not available and thus all known parametrix constructions for wave 
equations had to be restricted to domains V whose size is determined by the 
norm of the metric g in a given system of coordinates. Thus the Kirchoff-Sobolev 
representation would only hold for points p at maximal distance from E with 
dependent on the norm of the metric. As we shall explain below, such demand 
on the regularity of the metric would make the Kirchoff-Sobolev formula impossible 
to apply to realistic nonlinear situations, such as Einstein's field equations. 



Note that classically the construction of a Kirchoff-Sobolev parametrix could only be justified 
for points p such that J~{p) f] S belongs to a geodesically convex neighborhood I? of p. As we 
note below this requires uniform control for at lest two derivatives of the metric. 

^The method of [EM1],[EM2] was heavily dependent on the choice of a Cronstrorn gauge. 

''Defined as the image of the exponential map ; Tj,M — > M restricted to the largest convex 
subset of TpM where it is a diffeomorphism. 
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The importance of the classical condition becomes apparent upon examining the 
regularity of the null boundary N~{p) of the causal past J~{p). This set is ruled 
by past null geodesies 7(s) originating from p and terminating at the points 7(s*) 
beyond which one can find a time-like curve connecting p and 7(5) with s > s*, 
see [HE]. Regularity o{J\f~{p) breaks down precisely at the terminal points 7(5*). 
There are two reasons for the existence of a terminal point 7(5*). 

(1) 7(5*) is a conjugate point. 

(2) 7(5*) is a point of intersection of two different null geodesies. 

The existence of conjugate points is governed by the Jacobi equation for the Hessian 
D^u of the optical function u, 

Dl(D2u) + (D\f = R( • , L, ■ , L) 

with L = —g"^df3uda the null geodesic vectorfield along A/"" (p) and R the curvature 
tensor of g. This formula indicates that, at least as far as the conjugate points arc 
concerned, the terminal value of the afline parameter can be bounded below by 
an upper bound on sectional curvature which, in turn, can be controlled by a 
bound on the metric. 

Uniform bounds of of the curvature tensor R, or bounds for the metric g, are 

however not very useful in applications to nonlinear wave equations. For example 
in the classical local existence result for the Einstein vacuum equations [Br] , which 
is based on Kirchoff-Sobolev formula, the requirement is by itself worse^ when 
compared to the result in [HKM] based on the Sobolcv norm H^, s > 5/2. It 
is for this reason alone that the Kirchoff-Sobolev parametrix has been abandoned 
in all rigorous work on nonlinear wave equations in favor of energy estimates and 
Sobolcv inequalities. The main goal of our paper is to revive the Kirchoff-Sobolev 
parametrix by constructing it and showing that in the particular case of the Einstein 
vacuum equations, 

it is well-defined under much less stringent assumptions. For this task we rely in 
an essential way on the results in [Kl-Rol]-[Kl-Ro3] which show^ that the radius 
of conjugacy along J\f~{p), expressed relative to an afiine parameter of L, depends 
only on the size of the geodesic fiux of curvature® J^p along J\f~ [p) . These results are 
complemented by our recent work [Kl-Ro4] where we establish the remaining part 
of a lower bound on the radius of injectivity of Af~{p), i.e., control of intersecting 
null geodesies from p, expressed relative to a given time function. We achieve this 
by assuming, in addition to the above mentioned bound on the curvature flux, the 
existence of a coordinate system in V relative to which the metric g is pointwise 
close to the flat Minkowski metric. 



''^Additional losses of derivatives lead to a result in [Br]. 

^Properly speaking the results in [Kl-Rol]-[Kl-Ro3] do not consider the vertex p yet the meth- 
ods used in those papers can be shown to extend to cover the case of interest here. In fact, this 
forms the subject of the Q. Wang's thesis, Princeton University, 2006. 

^This is an appropriate integral of the tangential components of the curvature tensor along 
Af~{p), called curvature flux, which will be defined below. 
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2. Basic definitions and main formula 

2.1. Null cones. Consider a spacelike hypersurface S, a point p to its future S_|_ 
and J~{p) its causal past. We start by assuming the following local hyperbolicity 
condition for the pair (S,p): 

Al. All past causal curves initiating at points in a small neighborhood of J~{p) 
intersect E at precisely one point. 

Let Af~{p) be the null boundary of J~{p). In general J^~{p) is an achronal, 
Lipschitz hypersurface. It is ruled by the null geodesies^" from p, corresponding 
to all past null directions in the tangent space TpM. These null geodesies can 
be parametrized by fixing a future unit time-like vector Tp at p. Then, for every 
direction a; e S^, with denoting the standard sphere in M^, consider the null 
vector £^ in Tp(M), 

g(CTp) = l, (16) 

and associate to it the past null geodesic 7a; (s) with initial data 7a) (0) = p and 
7w(0) = ^u)- We can choose the parameter s in such a way so that L = 7a; (s) is 
geodesic. Thus, 

DlL = 0, g(L, L) = 0, and, at point p, g(L, Tp) = 1 (17) 

As mentioned in the introduction the null cone N~{p) is smooth as long as the 
exponential map {s,uj) Jlu{s) is a local diffeomorphism and no two geodesies, 
corresponding to different direction to €E>^, intersect. Thus for each w e either 
7a, (s) remains on the boundary of J'~ (p) for all positive values of s or there exists a 
a value s* (w) beyond which the points 7^, (s) are no longer on the boundary of J'~ (p) 
but rather in its interior, see [HE] . Thus TV" (p) is a smooth manifold at all points 
except the vertex p and the terminal points of its past null geodesic generators. 
Indeed, at a terminal point q there exists a null geodesic through q which fails to 
be in M~{p) past q. This implies that the tangent space Tq{M~ {p)) contains the 
past tangent direction of the null geodesic but not its opposite. This means that 
M~{p) must be singular at g. In what follows we shall denote by A/'~(p) the regular 
part of N~ (p) , that is the part with its terminal points removed. Clearly the null 
geodesic vectorfield L is well-defined and smooth on JV~{p). 



^'^Every point in A/" (p) \ {p} can be reached from p by a past null geodesic in Af (p). 
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The parameter s in the definition of is an affine parameter on J\f (p), i.e. 

L(s) = 1, sip) = 0. (18) 

Let 7 denote the degenerate metric induced by g on j{f~{p). Clearly 7(L, X) = 
for any X e TJ\f~{p). Let x denote the null second fundamental form of J\f~{p), 

x{X,Y) = s{T>xL,Y). (19) 

where X, Y are vector-fields tangent to Af~ (p) and D denote the covariant deriv- 
ative on (M, g). Clearly x is symmetric and x(L,X) = for any X G TJ\f~{p). 
This allows us to define trx as the trace of x relative to 7. 

Given a point q G N~ (p) \ {p}, we can define a null conjugate L to L such that, 

g(L, L) = -2, g(L, L) = g(L, L) = 0. (20) 

and further complement it by vectors (ei, 62) with the property that 

g(L,ea) = g(L, Ca) = 0, s{ea,eb) = Sab, a, 6=1,2. (21) 

The vectors (L,L, 61,62) can be locally extended to a neighborhood of a point 
q e 7V~(p) \ {p} to form a smooth local null frame. Relative to such a frame 

the only non-vanishing components of the null second fundamental form x S'^e 
Xab = g(De„L, Cb) = Xba- We Can introduce the other frame coefficients, 

X„^ = g(De„L,eb), Ca = ^g(D„L,L), = ig(e„, DlL) (22) 

Note that, in general, x^^b symmetric. 

Remark. A canonical way to define a null geodesic conjugate is to take L the 

unique null vcctorficld orthogonal to the level surfaces Ss defined by the affine 
parameter s. We refer to the corresponding null pair as a null geodesic pair. We 
can also choose ei, 62 to be tangent to Sg- Note that in that case x is symmetric. 
We also note that in a neighborhood of p where (p) coincides with its regular 
part Af~ {p) the geodesic null frame defined above is smooth away from the point 
P- 

For the purpose of constructing our Kirchoff-Sobolev parametrix we shall make, in 
addition to Al the following assumption. 

A2. We assume that J\f~ (q) coincides with M~ {q) past the space-like hypersurface 
S for any point q in a neighborhood of p. 

2.2. Optical function. To make sense of our Kirchoff-Sobolev formula we need to 
define an optical function^ ^ u, in a neighborhood of M~{p), such that it vanishes 
identically on JV^{p). We define u uniquely relative to the time-like vector Tp as 
follows: 

Let e > a small number and : (1 — e, 1 -|- e) ^ M denote the timelike geodesic 
from p such that Te{l) = p and r^(l) = Tp. From every point q of Fe let M~{q) 



i.e. a function which verifies (11) 
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be the boundary of the past set of q. In view of assumption A2 for ah sufficiently 
small e > 0, J^~{q) coincides with its regular part J^~{q) to the future S+ of S. 

We now set u to be the function, constant on each Af^{q), such that for q = T^{t), 

This defines a smooth function u which vanishes on J\f~ (p) and verifies the eikonal 

equation (11) 

in a neighborhood of JV^(p) D S+. Observe that the null geodesic vcctorfield 
L = ^"'^dfjuda extends the vectorfield in (17) to V. It verifies the normalization 
condition, 

g(L,Tj,) = Tp(u) = l, 

at all points of T)^. We can thus extend the definition (19) of the null second 
fundamental form x and its trace trx at every point in D^. 

We can introduce local coordinates around any point in r € by considering the 
unique null geodesic ^uj,q, with € S^, which initiates at 9 € Fe and passes through 
r at value s of its affine parameter. Denoting by u the value corresponding to the 
null cone J\~{q) we see that r is determined by the coordinates u, s and w e 

2.3. Dirac measure on N^{p). Given our smooth optical function u, defined 
in the neighborhood T), of Af^ (p) H S+, and a distribution /i on the real line 
R, supported at the origin, we can define the pull-back distribution u* (fx) — fx o 
It on C M in the usual sense of distribution theory. In the particular case 
when /i is either the Dirac measure So or its derivatives S'q,6q, . . . , we denote the 
corresponding distributions on M by S{u), 6'{u), S"{u),... . We can thus make 
sense of calculations such as, 

D„(5(u) = S'{u)'DaU, DcD/3(<5(u)) = S"{u)'DaU'Df3U + (5'(u)D„D/3U 

Clearly S{u), 5'{u),. . . are supported on 7V~(p) n T>^. We can use the definition of 
S{u) to define the integral along N~{p) of any continuous function / supported in 
as follows. 

Definition. Given a continuous function / supported in we define its integral 

onM-(p) by, 

/ f=<5{u)J> (23) 

Jjir-{p) 

Proposition 2.4. The definition (23) depends only on the restriction of f to N'{p) 
and the normalization condition (16) used in the definition of the null geodesic 
generator L. 

Proof : We may assume without loss of generality that / is supported in the 
domain V^, which can be parametrized by the coordinates u, s and w G as 
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described above. We can the easily calculate, according to the definition of d{u) 
and coarea formula, 

<Siu),f>= f f fiO,s,Lu)d.sda, 
Jo Js^ 

where dus denotes the area element on the 2- surfaces Sg of constant s. ■ 

2.5. KirchofF-Sobolev parametrix. Consider Jp to be a fixed k-tensor at p and 
let A be the unique k -tensor-field defined along J\f~ (p) which verifies the linear 
transport equation, 

DLA+itrxA-0, {sA){p)=3p (24) 

with s the afiine parameter (18). The tensor-field A can be extended smoothly^^ 
to a small neighborhood of M^{p). We can now define the distribution, or current, 
in E+, 

< A5{u),¥ >=< S{u), g(A, F) > (25) 

for an arbitrary, smooth, fc-tensor- field F supported in S+. Here g(A,F) denotes 
the full contraction of the k -tensor-fields A and F with respect to the space-time 
metric g. Observe that the current A5{u) depends only on the choice of Tp and 
Jp and not on the particular extensions of u and A. 

In what follows we identify the space of /e-tensors at p and its dual with the help 
of the metric g. 

Definition. We call ICp = K^pj^, a k-tensor- field distribution with values in 
the space of k-tensors at p, defined by the formula ^p = Ad{u), with A defined 

by (24), the retarded Kirchoff-Sobolev parametrix at the point p, corresponding to 
Jp. If "il is a solution of the equation □gVE' = F, with F supported in S+, we denote 
by 'S'F,jp(p) the k-tensor at p defined by the integral, 

*F,j,(p)=</C-j^,F>= / g(A,F). (26) 



In the case of the scalar wave equation Ogijj = f wc can choose A to bo the scalar 
solution of (24) with initial data {sA){p) = 1. In that case we have JCp = A5{u) 
and 

^f{p)=<K-,f>= [ Af. 

In the particular case of Minkowski space we can easily identify A = Ap{q) with 
the term |a; — y\~^ where q = {s,y) e M~{p) and p = {t,x). 

"'^^We can in fact extended it canonically by solving the same transport equation along J\f^(q), 
with q € Fe and initial data sA{q) = Jq where Jq is an arbitrary smooth tensor-field coinciding 
with Jp at p = q and s the afine parameter along Af~{q). Note that, so defined, the tensor- field 
A is smooth away from the axis T^. 
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2.6. Time foliation near vertex. Returning to the construction of u in subsec- 
tion (2.2) we observe that the parameter t along the geodesic Fe can be extended 
to a local, equidistant^^ time foliation T,t, t £ [1 — e, 1 + e] which covers a whole 
neighborhood of the point p, such that p G Ei. Indeed, starting with a fixed space- 
like hypersurface Si through p, orthogonal to the future unit timelike vectorfield 
Tp, we can define this geodesic foliation using the timelike geodesies normal to Si. 
In particular, for all t € [1 — e, 1] , if we denote by fie the set 

= {J- (p) n S+) \ U,e[i_,,i] St (27) 

then its boundary is given by 

afie =A/'-(p)U£)i_eU£) 

where M~{p) is the portion of M~{p) to the future of S and the past of Si_e, 
Di-e = J-{p) n Si_e and D = J-{p) n S. 

Let T = Dt denote the future, unit normal to the foliation S^, defined in a neigh- 
borhood of p. We define the null lapse function (p and the second fundamental form 

k associated to S^: 

<^-i =T(u) = g(L,T), fc(x,y) =g(DxT,r), VX,reTSt. (28) 

Clearly (p{p) = 1. Since T is a locally smooth vectorfield, k is a smooth symmetric 
2-tensor. In particular, 

||fc||L=o <C 

for some constant C. Similarly, since u is a smooth optical function and (p{p) — 1, 
the lapse <^ is a smooth bounded function in a neighborhood of p. in particular, 

|^((7)-1| -»0, q^p. (29) 
We now recall the Raychaudhuri equation satisfied by trx along X^(p), 

^(trx) + ^(trx)' = -\X? - Ric(L, L). (30) 
with s the afine parameter of L and x the traceless part of x- 

The behavior of the function trx at the vertex p is determined by the conditions 

(strx)(p) = 2, X(P) = 0. (31) 
Integrating the Raychaudhuri equation one can easily deduce that, 

|trx(9)-^H0, q^p. (32) 

Consider the time function t restricted to M~{p). Then 

-=L(i)=g(L,T)=<p-i (33) 
The area \St{p)\ of the 2-d surfaces St{p) = St n A/'~(p) obeys the equation 



d f 

-T:\St(p)\ = / iftrxda^. 



^^With the lapse function of the fohation identically one. 
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This and the behavior of tr^ and (f> near p {t{p) = 1) imply that 

\St{p)\=Mt-if + Oi\t-l\') (34) 
On the other hand from (33) and (29), t — 1 = s + o(s), which implies that 

\Stip)\=Ans^ + o{s^) (35) 

We shall also make use of the following simple variation of proposition 2.4. 

Proposition 2.7. Let t be a regular time function defined onM~{p) with t{p) = 1 
and equal to <1 on YioC\N~ {p), where Sq is an arbitrary spacelike hypersurface on 
J\~{p) n S+. Assume that <^ = || < 0. Then, for every test function f , compactly 
supported in J^{T,o) 



-fd. 



<S{u),^>= I j fifdtdat (36) 
where St denotes the level surfaces of t and dat the corresponding area element. 



Proof : The result follows easily by first extending t and u to a neighborhood D 
of Af~ (p) n S+ and the applying the coarea formula as above. ■ 



2.8. Statement of the result. We consider a space-like hypersurface E c M and 
a point p e S+ = J7'+(S) such that the assumptions A1-A2 are satisfied. 

Theorem 2.9. Let 6e a solution of the equation □g^' = F with F a k-tensor-field 
supported in Then for any k-tensor 3p at p, 

*(p) = *F,j,(p)+ / g(^,*), (37) 

where 

*F,j,(p)=</C-j F>= / g(A,F) 

and A verifies (24). The smooth error term £ depends only on Jp, the geometry of 
the truncated null cone J\f~ (p) D S+ C M, and the ambient spacetime curvature R 
restricted to N~ (p) . 

In the particular case of a scalar wave equation O^tp = f, A and £ are scalar 
functions on N~ (p) and, 

V'(P) = V'/(P) + / ^V', 



The precise expression of the error term £ will be given in Theorem 3.11 
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3. Derivation of Kirchoff-Sobolev formula 

3.1. Covariant derivatives of space-time tensors. As is well known there is no 
canonical way to define a restriction of the space-time covariant derivative D to a 
null hypersurface. This is due to the absence of a canonical projection of a tangent 
space TqM, q € Af~{p), onto the tangent space Tq{Af~ (p)). This projection can 
be fixed, however, by a choice of a null conjugate L, i.e. a null vector such that 
g(L,L) = —2. With this choice we define an induced covariant derivative ^—^D on 

^^DxY = T>xY + ^xiX, Y%, MX, Y e TN-{p) 
For example, if we choose X, Y to be the elements {ea)a=i,2 of a null frame (L, L, ei, 62 

We now make sense of covariant derivatives of space-time tensors along J\f~{p). 
We start by defining a covariant derivative D of a space-time 1-form defined 
on j{f~{p). Thus wc view ^ as a section of the vector bundle "'^'^ T*M over J\f~{p), 
endowed with the induced covariant derivative ^D. We interpret the covariant 
derivative T>A of A along J\f~ (p) as a 1-form on TV" (p) with values in T*M. Thus, 
for every vectorfield X G TN~{p) and any vectorfield Z in TM, 

tiA{X;Z) = t>xA{Z) ■.= X{A{Z)) -A{T>xZ) 

We also write, 
" 2 

Wc define D A, an Af~{p) 2-tcnsor of second covariant derivatives of A along 
M~{p) with values in T*]V[, by the formula, 

D^A(X, Y-Z) = { Dx DA)(y; Z) = X{ t)A{Y; Z)) - tiA{^^^DxY; Z) - hA{Y- D 
or simply, 

D%(X,r) = (Dx(f)yA)V - (D(i)z5^r^)^ 

These definitions can be easily extended to higher covariant derivatives along M~ (p) 
and to higher order tensors A. 

3.2. Kirchoff-Sobolev current. Consider the current fC~j^ = AS{u) defined in 
(25). Recall that Jp is an arbitrary fc-tensor at p and A is a fc-tensor- field verifying 
the transport equation, along M~{p), 

DlA + ^Atrx = 0, sA{s)\s=o = Jp- (38) 

Also L = g^"dvudfj,, gi^^dvudfjiU = and L(s) = 1. Let L be an arbitrary local 
null conjugate to L, i.e. JAu) = g(L,L) = —2. Calculating relative to an arbitrary 



with the covciriant derivative denoted by D 
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null frame wc easily check that DgU — tvx- Formally we thus have, 
□g(A5(w)) = □gA5(w) + (g'"^D^AD,.M + Aagu)5'{u) + A{g'"'d^ud^,u)d" {u) 
= DgA6{u) + (-L(m)DlA + Aagu)d'{u) 



□gA(5(u) + 2(DlA + -Atrx)^'(u) 



Hence, 



□g(A,5(«)) = □gA5(u) + 2(DlA + ^Atrx)S'{u). (39) 



Observe that the above calculation does not depend on the choice of L. 
Since 

DL(sA) = -i(trx-^)sA (40) 

we have, in view of (32), that along N'~{p), 

|sA(g)-Jp|^0, s^O. (41) 

We shall next apply /C~j = A5{u) to the equation Dtp = F in the sense of 
distributions, 

/ g(A^(n),ng*)= / g{a{AS{u)),F) (42) 

where S+ = J7'"'"(S) is the future of the initial hypersurface S. We assume that 

has zero data on E and that F is supported in E"*". Our next goal is to integrate 
by parts on the left hand side of (42). We first decompose DgA, for an arbitrary 
tensor-field A, relative to our null frame (20) - (21). For simplicity we assume that 
A^ is a one tensor, the general case can be treated in the same manner. We recall 
the definition of the Ricci coefficients (22), 

X„t, = g(DealL,eb), Ca = ^g(DaL,L), = ^g(e„, DlL) 

and also introduce, 

a; = -Jg(DLL,L) (43) 

which is well defined in a neighborhood of N'~ (p), see subsection 2.2. Using also 
the notation in subsection (3.1) we derive: 

Now, DfeA^ = t>bA^ and, since JDaCb = ^-^DaCb + ^XabL., 

^lb\ = ea(DfcA^) -DD^e.A^ 

= ea( Db^^) - D(L)£,^e,^/^ - ■^Xab'Di^Ai^ 
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Hence, denoting AA^ = S"''' t)l^A^, 
On the other hand, 

= DlDl - 2Ca t>aA^ + 2uj DlA^ + R;^ ll 

Henceforth, 

□gA^ = -DlDlA^ + AA^ + Ca ■ t>aA^ (44) 
- uty^Ai, - ^trxDL^^ - ^R^^ LL^A 
Remark 3.3. In the case when A is a scalar formula (44) becomes, simply, 

□gA^ = -DlDlA + AA + Ca- t>aA - Lo Dl^ - ^trxDj^A (45) 



3.4. Integration by parts. In view of (39) we have, 

□g(A(5(u)) = □gA5(M) + (2DL^ + trxA)5'(u) (46) 
where for u = 0, 

DlA + ^Atrx = 0. 

According to section 2.6 we have defined a time foliation Sj, t G [1 — e, 1 + e], in 
a neighborhood of the vertex p with p g Si such that the boundary of the set 

= J^ip) \ Utg[i_e,i] is given by, 

dVL, =jif-{p)UDi^,UD. 

Here Af~{p) is the portion of M~ (p) to the future of Si_<; , -Di_£ = J~{p) n 
and D = J~{p) n S. As before we denote by T the future unit normal to the 
surfaces Sq = S and Ei_e. Note that Ei_e needs only be defined locally, for small 
e > 0. Note also that, for e = 0, T coincides with Tp as defined in section 2.1. 
Clearly, 



/ g(A%),ng*) = lim / g(A5(,i),ng*) 



where 17 = J^{p) H S+. Due to the presence of 5{u) and the fact that ij} is 
supported in S+ we may assume in what follows that all functions we deal with 
in the calculation below are supported in the set Vl = J^{p) H J+(E). Thus the 
boundary of the intersection of their supports with is included in Ei_e. 

Lemma 3.5. Let F, G he two tensor-fields of the same rank and F is a distribution 
supported in fi. Then 

t=l-e 



I g{F,D^G)= [ g{n^F,G)- f (g(F,DTG)-g(G,DTF)) 

JQe -'fie •'Dt 



t=0 



where Dq = D. 
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Proof : Indeed, 

g(F, DgG) - g(ngF, G) = D"g(F, D,G) - D"g(D,F, G) 

Thus, 

/ (g(F,ngG)-g(ngi^,G)) = / D"(g(F,D«G) - g(D„i^,G)) 



Dt 



T"(g(F,D„G)-g(D„F,G)) 



l-e 



We now write, 

g(A(5(u),agVi/) 



g(ng(A(5(u)),«') 



- / g(A5(u),DT*) + / g(DT(A5(u)),*) 
= / g(ng(A(5H),M/)+/, + Je, 



where 1^ and Jt denote the boundary terms on -Di-e. The term corresponding to 
Da vanishes due to the zero data assumption for 

Proposition 3.6. We have 

le^O, Je ^ -47rg(\E'(p),Jp) as e ^ 0. 

Thus, 

/g(A5(«),ng*) = hm / g(Dg(A5(w)),*)-47rg(*(p),Jp) (47) 



Proof : We analyze the boundary terms Ie,Je- Clearly, 

Ie = - [ g(A^(w),DT*) = - / g(A,DT*)^da^ 

where, see (28), (p = |Dtw|^^ is the null lapse and da^ is the area element of 
the 2-surface Si-e{p) = ■Af~{p) n Di-e- Recall that according to (34) the area 
\S^.,{p)\<e\ 



Now, 



|/,| < ||^|Uo.( / \A\'da^)'^\ f \T>T^'da. 

< ||^|U^||DT*|U^||A|U.(s,_.(rt)|5i_,(p)|i/2 

< e ||95||L-||DT*||L-||A||i2(5^_^(p)) 

Recalling (41) and (34) we easily see that ||A||i2(Sj_^) is bounded as e — > 0. Thus, 
for a smooth tensor-field ^, we clearly have. 
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We now consider the second boundary term, 

Je = [ g(DT(A(5(«)),vl/) 

= / 5(u)g(DTA,*)+ / S'{u)r>TUg{A,^) 
= [ %)g(DTA,*)+ / 5'(«)(p-ig(A,*) 

= J} + J! 

UN — ipL + T denotes the unit nornial"'^''' to S'i_e = Af~{p) D Si_e in then 
Djv(5(m) = (5'(M)DjvM and D^ru = Dtw = </?^^- Hence, 

J2= / <5'(w)<^-ig(A,f) = / D^5(w)g(A,*) 
We next record the following integration by parts formulae. 

Lemma 3.7. Let X be a vectorfield tmigent to the hyperplane St and let f,g be 
two scalar functions on Ej. Denote by V the covariant derivative restricted to S^. 
Then, 

[ fX{g) = - [ {X{f) + div Xf)g. (48) 

In particular, 

I fN{g) = ~ [ {N{f) + tref)g, (49) 
where tr6 is the mean curvature of the 2-d surfaces St,u = {u = const} fl St. 

Proof : Formula (48) is standard. To prove (49) observe div TV = g(VArA'', N) + 
J2a g(^a-^> ^a) = tr0 where is the second fundamental form of the surfaces St,u C 
St. ' ■ 

Using the lemma we infer that, 

J! = - I <5(u)(JVg(A,*)+tr^g(A,*)) 

= - / ,5(«)(g(D^A,*)+tr0g(A,*)) - / %)g(A,D^*) 

Now, proceeding as for 7^, it is easy to check that Jj, (5(u) g(A, Djv^) ^ as 
e ^ 0. Hence, 

limj2 = -lim / (5(u)(g(DjvA,*)+tr6'g(A,*)) 

Or, 

limJe = -lim / ^(M)(g(Djv-TA,*)+tr^g(A,*)) 



^^A. priori, the vectorfield N is defined only on J\f (p) n Si-^, it can however be extended 
locally as a unit normal to the foliation of 2-d surfaces {u = const} fl Si_e. 
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Now observe that L = if~^{N - T). Hence, DlA = (^"^Djv-tA. Since DlA + 
^Atrx = we infer, 

Djv-tA = vDlA = -^ipAtvx- 

Hence, 

lini Je = - lini j 6{u){ - ^(ptYX + tr6')g(A, *) 

On the other hand 9ab = gC^aN, Cb) = g(Da(<^L + T), eft) = ^pXab + Kb- Therefore, 
tr^ = v'trx + S"'^kab and we deduce, 

hm Je = -I lim / ^tvxS{u) g( A, - lim / {S''''kab)6{u) g( A, -J-) 

It is easy to see that the second term of the right hand side converges to zero for 
e ^ 0. Indeed, 

\[ (5»%fc)5(w)g(A,M/)| = \f (5-''kab)s{A,^)^da^\ 

< |5i_e(p)|^||fc|U~||vI/|U»||<pU»||A|U.(5,_.(rt) 

Therefore, 

hniJe = -l-lim [ (/j^ trxg(A, da^ 

It is easy to check that, 

I Urn / ^j.^ g (* - da^ = 

Therefore, 



hniJe = -ihm/ (^^trx g (A, da^ 



Or, since sups^_^(p) - 1\ ^ 0, and sup5^_^(p) |trx - f | ^ as e ^ 0, we infer 
that 

hmJj = -e~^Hni/ g (sA)) = -47r hm g (^'(p), (rA)(r)) 



hm / 



Thus, using the initial condition Hmg^o sA{s) = Jp, we obtain 

hniJe = -47rg(*(p),Jp) 

e— >0 



We now analyze the term g (□g(A(5(u)), on the right hand side of (47). In 
view of (46) we have, 

/ g{a^{AS{u)),^)= [ 5(w)g((ngA),*)+ / ,5'(«)g((2DLA + trxA),'J'). 
Jq. Jq, Jn, 
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Given the normalization L(u) = —2 we have S'{u) = — iDLi5(u). Integrating by 
parts we obtain 

/ y(7/)g((2DLA + trxA),*) = i / (5(u)g(DL(2DLA + trxA),*) 

+ i ^ 5(u) g ((2DlA + trxA), (Dl* + D«L„*)) 

+ / 5{u) g(((2DLA + trxA), *) g(L, T) Ijlj"' 

Recall that 2DlA + trxA = on the surface it = 0. Therefore the last two terms 

vanish and wc derive, 

/ 5'(u)g((2DLA + trxA),M/)= / g f Dl(DlA + ^trxA), * 
Therefore, 

/ g(ng(A^(w)),*) = / ,5(ix)g((ngA + DL(DLA+itrxA)),vl/') 
We now recall (44), 

□ gA = -DlDlA + AA + Co • DaA 



- wDlA- ^trxDLA- ^R(-,A,L,L) 



Therefore, 



□gA + Dl(DlA + itrxA) = AA + Ca f>aA - ^trxDL A 

- w Dl A + 1 (DLtrx) A - iR(., A, L, L) 



while since DlA + ^trxA = on TV (p). 



-itrxDLA - ijDlA + i (DLtrx) A = ^(DLtrx + ^trxtrx + 2wtrx)A 
Hence, we have proved the following. 

Proposition 3.8. In the case of a one form A verifying (38), 
1 

2' 



□gA + DL(DLA+-ir^A) = AA + CoD„A 



+ \{^htrX + \trxtrx + 2a;trx)A - A, L, L), 

where 

AA = AA + ^ir^DLA = e„(D„A) - g(D„e„, 6^)05^ 

zs a Laplace- Beltrami type operator which coincides with the standard surface Laplace- 
Beltrami operator in the case when the frame {ea}a=i,2 spans a tangent space of a 
2- dimensional surface^^. 



in the case of the geodesic foliation. 
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In the scalar case we have instead, see remark 3.3, 

□gA + DL(DLA+iirxA) = AA + CaDaA 

+ \{^hP'nc + ^trxtrx + 2u}trx)A, 

Using the above proposition we infer that, 

/ g(Dg(AJ(«)),*) = / 5(«)g((AA + CaD„A),*) 

If 1 

+ - / 6{u) (DLtrx + ^trxtrx + Swtrx) g(A, *) 

+ i / 5(w)R(*,A,L,L) 
We now make use of the following, 

Proposition 3.9. Introduce the mass aspect function as in (13.1.10b) o/[C-K], 

fj, = 'DiJrx+^trxtrx + 2u;trx, (50) 
The following formula holds true relative to the standard geodesic foliation on 

f, = 2divC-X-X+ 2|Cl' + Rll + ^R(L, L, L, L) (51) 

Proof: See [C-K]. ■ 

Remark 3.10. Note that according to (51) the mass aspect function n depends only 
on the null hypersurface A/"" (p) and the ambient curvature R. 

We have therefore proved the following precise version of theorem 2.9, 



DlA + -Atrx = 0, sA{p) = 3p on u = 



Theorem 3.11. Let A be a vectorfield verifying, 

1 

2' 

where 3p is a fixed vector at p. Then solution ^ of on inhornogeneous vector equa- 
tion Dg^* = F, in a globally hyperbolic spacetime (M,g) satisfying A1,A2, with 
zero initial data on a Cauchy hypersurface S can be represented by the following 
formula at point p with fl = J~{p) n J7'+(S), 

47rg(f(p),Jp) = - / <5(«)g(A,F) (52) 

-If 5(m)R(v1/,A,L,L)+ / %)g((AA + CaDaA),*) 
^ Jn Jn 

+ 11 S{u)ng{A,^) 



n 



2 

where, 

R(*, A,L,L) = R„/3^5LTL^*«A'^, 
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with 'Rafj-yS the components of the curvature tensor R relative to an arbitrary frame. 

Remark 3.12. Theorem 3.11 implies that the error term £ in (37) has the foUowing 
representation 

£a = -iRaA^^A^L^L* + (AA + CaD„A)„ + |a„ 
in the case of a vectorial wave equation. For the scalar wave equation 

£ = (AA + CaD.A) + |a 

Remark 3.13. Formula (52) can easily be generalized to higher order tensor wave 
equations. Indeed if both ^ and F are tensor-fields of order k then Jp, A and £ 
are also of order k and, 

g(f,*) = g((AA + CaD„A),*) + |g(A,*) (53) 

+ R(-,-,L,L)#*#A 

where the last term denotes a scalar contraction of R(-, -jL ,L) with and A. 

4. Wave equation for sections of vector bundles and applications to 

THE Yang-Mills equations 

Now let V be a vector bundle over (M, g) with a positive definite scalar product 
<, > and a compatible connection A. We may assume that V is a vector bundle 
associated to a principal bundle P so that V = P xq E with G a compact Lie 
group and a vector space E. Let G denote the Lie algebra of G. The connection A 
is a ^ valued 1-form on V, which, locally can be viewed as a ^ valued 1-form on 
M. 

We define the gauge wave operator for sections ^ : M — > V 

where = + [A^, ■ ] denotes the gauge covariant derivative. We denote by A 
the curvature of the connection, i.e. the G valued 2-form on M, 

Aa/3 = daX/3 - dpXa + [Aa, A/?]. 

As before we construct a Kirchoff-Sobolev parametrix /C^ for □g'^ by defining 

K~ = /C~T = AS(u), where A is a section of V which verifies the covariant 
transport equation^ ^ 

DlA + ^trxA = (54) 

with initial data (sA)|,5=o = Jp and Jp is a fixed clement of the fiber Vp. As before 
we assume that (M, g) is globally hyperbolic and satisfies Al, A2. We also assume 
that u is a solution of the cikonal equation g^^^daudisu = with u, vanishing on 



Note that here transport of A along integral curves of L is modulated by the action of the 
gauge potential A. 
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the boundary M [p) of the past of p in M. Repeating our calculations of section 
3, we obtain the following analog^^ of Theorem 3.11. 

Theorem 4.1. Let A he a section of a vector bundle V over (M,g) verifying, 
DlA + ]^Atrx = 0, sA{p) = Jp on J\f-{p) 

where 3p is a fixed element of\p. The solution \E' of the inhomogeneous gauge 
equation Dg^^ = F, with zero initial data on a Cauchy hypersurface S can be 
represented by the following formula at point p with ft = J~{p) n J7'"'"(S), 

47r<*(p),Jp> = -[6{u)<A,F> (55) 
Jn 

- \ [ 5{u) < [All,A],vI/ > + / 5{u) ((A(^)A + Ca2?i^^A),vI/^ 

+ \ 5{u) II <A,^> 
^ Jn 

In particular, in Minkowski space, for general initial data, we have the following 
representation: 

47r < *(p), Jp > = - f S{u) <A,F>-\ I 5{u) < [All, A],* > (56) 
Jn ^ Jn 

+ I S{u) < A(^) A, * > + / (< AS{u), Dt* > - < T>t{AS{u)), * >) 
Jn JT, 

The last term represents contribution of the initial data on S. 

Remark 4.2. The formula (55) can be naturally extended to consider sections 
of the bimdlc TM ... (8) TM V. In this case terms of the form 6{u) < 
R(-, ■, L, L)A, \l/ >, where R is the Riemann curvature tensor of g, need to be 
added. The corresponding extension of (56) does not therefore introduce any ad- 
ditional terms. 



4.3. Yang-Mills equations. We now assume that A is a Yang-Mills connection 
on a 4-dimensional Lorentzian manifold (M, g), i.e., it verifies the equations 

2?«A„^ = 0. (57) 

The Yang-Mills equations are hyperbolic in nature and admit a Cauchy formulation, 
in which the connection A is prescribed on a Cauchy hypersurface^^ S and then 
extended as a solution of the problem (57). The uniqueness and global existence 
for the Yang-Mills equations with smooth initial data in 4-dimcnsional Minkowski 
space-time was established by Earcllcy-Moncricf, [EMI], [EAI2]. Tliis result was 
later extended to the Yang-Mills equations on a smooth 4-dimensional globally 
hyperbolic Lorentzian space-time by Chruscicl-Shatah, [CS]. A different proof in 
Minkowski space, allowing for initial data with only finite energy , was given by 
Klainerman-Machedon, [KM] . 

"'^^Thc only new term in the formula (55) below is due to the commutator (D^^D^"^' — 
D^'^'d^^')A = [AiL, A] in deriving formula (44) 

^^This requires space-time M to be globally hyperbolic. 
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All of the above approaches were manifestly non-covariant; as they required a choice 
of a gauge condition for the connection A. The approach of Eardley-Moncrief was 
based on the fundamental solution for a scalar wave equation in Minkowski space 
(Kirchoff formula) and made use of Cronstrom gauges: For any point p the connec- 
tion A can be chosen to satisfy the condition 

The work of Chrusicel-Shatah relied on the Friedlander's representation of the fun- 
damental solution of a scalar wave equation in a curved space-time and a local 
analog of the Cronstrom gauge. Finally, Klainerman-Machedon's proof was based 
on a Fourier representation of the fundamental solution of a scalar wave equation 
in Minkowski space, bilinear estimates and the use of the Coulomb gauge: 

au, 

Below we present a new simple gauge independent proof of the global existence and 
uniqueness result for the 3 -|- 1-dimensional Yang-Mills equations. The main new 
ingredient is the use of a gauge covariant first order Kirchoff-Sobolev parametrix 
described in Theorem 4.1. 



Differentiating the Bianchi identities ©[^A^^.] = and using the equations we infer 
that the curvature A is a solution of a covariant gauge wave equation 

□^^)A„^ = 2[A"„, A,,^] + 2R^„7/3A"'' + Ra.A/ + R^j^A"^ 

For simplicity we consider the problem in Minkowski space, although our results can 
easily be extended to the general case of a globally hyperbolic smooth Lorentzian 
manifold, as in [CS]. The equations then simplify, 

□ (^)A„^, = 2[A-„ , A^p] (58) 

Recall that the curvature A can be decomposed into its electric and magnetic parts 
Ei = Aoi and Hi = *Aoi. We also recall that the total energy 




is conserved. Moreover, adapting the energy identity to to the past J'~ (p) of a point 
p G we also get a bound on the flux of energy along M~{p). More precisely we 
derive, 

where the backward null energy flux is defined with the help of a null frame 
(L,L, ea) centered at p. Without loss of generahty we may assume that p = {t,0) 
and denote r = \y\. Then L = dr—dt,^ = —dt — dr and Ba is a frame on a standard 
sphere §r- With these notations 

/ (iALLi'+Ei^L«n 

As in the original approach of Eardley-Moncrief the key element of the proof of 
global existence is a pointwise bound on curvature A. Once this bound is established 

the remaining steps concerning existence, propagation of rcgiilarity and uniqueness 
are very standard and will be omitted. The precise statement concerning an L°° 
bound on A is as follows: 
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Lemma 4.4. There exists t* > dependent only on So such that for any point 
p = {t, 0) we have 

\Mp)\<Ct-r„ 

where the constant Ct^r, depends only on the solution K on a hypersurface St_T-^ . 

Remark 4.5. Iterations of Lemma 4.4 leads to a pointwise bound on the curvature 
A in terms of the initial data. 



Proof : Wc fix r* > 0, whose is to be determined later, and apply the represen- 
tation formula (56) in the domain fl = J^{p) H {Y,t-T^), 

47r<A(p),J> = -2 / 5{u) <A,[A,A]>-i / 5{u) <[All,A],A> 
Jn ^ Jn 

+ / 5{u) < A(^)A,A> 
Jn 

(< A5{u), DtA > - < Dt(A5(u)), a >) (59) 



Here J is an arbitrary Q valued anti-symmetric 2-tensor on M^+^, A is a ^ valued 
2-form on M^+^ verifying the equation^" 

2?LA + r-iA = 0, (rA)|^=o=J 

and <, > denotes a positive definite scalar product on A?(M?^^)®Q. The last term 
in (59) depends only on the solution A on T,t-T, and therefore is consistent with 
the claim of Lemma 4.4. We now observe that for a, 6 € A^(M^+^) ^ we have 

\ < a,b > \ <\a\ \h\, 

where \a\ denotes the absolute value of an element in A^(R^) ® Q relative to the 
positive definite scalar product^^. In what follows all the norms will be understood 

to involve the absolute value | • | on A^(]R'') (g) Q. We denote by N~{p) the null 
boundary of O to the future of . Then 



J{u) < A, [A, A] > I < ||rA||^„(^-(^)Jr-i[A,A]||^,(^-(p)), 
I / 6{u) < [All, A], A > | < \\rM\L^(M~(p))\\r~'^-Li,\\Lii^-(p))\\^^\\L^iM-(p))' 



1/ 

Jn 



J{u) < aWa, A > I < ||A||^^(^- (^))||AWa||^,(^- 
It is easy to see, see e.g. [EM2] that 

|[A,A]|<|A| (iAllI + EI^l-I) 



a=l 



and therefore 



Similarly, 



r 



"'[A, A]||^,(^-(^)) < r.' ||A||^„(^-(^„ {J^-y 



\r 'ALL|lii(A^-(p)) (^p )■ 



20Recall that L = dr - dt, = dr - dt + [Al, •] and r = \y\. 
■^^JR^ here stands for the Euclidean 4-dimensional spa<;e. 
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To prove Lemma 4.4 it would be sufficient to show that 

II^A|l^~(Ar- W) ^ < |J| (rl {:F-y+nJ^-) 

(60) 

and then choose r, << (1 + £o)~^ as < Eq. 

To prove (60) we introduce a new Q valued 2-form B = rA so that 

I?lB = 0, B|,=o=J- 

Considering the components of the 2-form B it suffices to assume that B is a 5 
valued function on 1R^+^, in fact on J\f~{p). 

Commuting^^ the transport equation PlB — with r^A^^^ we obtain 

Pl(t-'A(^)B) =r2[AL",Vi^)B] Wi^\A^,B] 
We also have the equation 

PL(rVi^)B) = r[ALa,B] 
We combine these equations into the system: 

PlB = 0, (61) 
pL(rVi^^B) =r[ALa,B], 

I?L(r2A(^)B) = 2r2 [Al", V^^'B] + [Vi-^)AL'', B]. 
The first equation immediately implies that 

sup |B| < |J|, 

as the covariant derivative D,,, = da + [A.-] is compatible with a scalar product on 
G. We infer from the second equation that 

||V(^)B|U=(5,) < |J| fil \^^a?daA dp 

where dcjs is the are element of a 2-dimensional sphere Sp of radius p. 

To treat the last equation in (61) we need to worry about the term [v1'''^Al,B] 
which contains derivatives of A. Recall that the flux only allows us to estimate 
the tangential components of A and none if its derivatives. We get around this 
difficulty by expressing the Yang -Mills equations D^A^/j = relative to the null 
frame L,L, ei, 62- In particular, D^Alo + 'P— All = 0. This in turn implies that 

V^-^^Al" = JplAll + -All. 

Thus 

pL(r^A(^)B - ir^[ALL,B]) = 2r^ [K^MaM 

^^Recall that V*^"**^ is a gauge covariant derivative acting on sections Sr ^ P y.AdS and A''*') 
is the corresponding gauge Laplace-Beltrami operator on a standard sphere Sr of radius r. 
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Therefore, 

||AWB|U.(s,) < |J| ^i||ALL|U.(5.) [ \^i.a\'da)j dp 

Integrating with respect to r we obtain 

and the result follows. 



0=1,2 / 



5. Applications to General Relativity 

In this section we specialize our results to Einstein vacuum space-times (M, g): 

R„0-iRga;3-O (62) 

Equations (62) combined with the Bianchi identities imply that the Riemann cur- 
vature tensor Ha^i^u of an Einstein vacuum metric g satisfies a covariant wave 
equation 

where the quadratic term R ★ R is obtained by a contraction of the curvature 

tensor HaPfiu with itself. 

Theorem 5.1. Let p be a point to the future of a space-like hypersurface E in 
an Einstein vacuum space-time (M,g). We assume that assumptions A1,A2 are 
verified at p. Let A be a 4-tensor verifying, 

DlA -I- ^Atrx = 0, sA{p) = 3p on u = 

where 3p is a fixed 4-tensor at p. Then the curvature tensor Hap^v of g can be 
represented by the following formula at point p with O = J~{p) n J^{Y?), 

47rg(R(p),Jp) = - / (5(w)g(A,R*R) (63) 
- \ j ,5(«)R(-,-,L,L)#R#A+ / <5(n)g((AA + CaD„A),R) 

+ \ l 5(«)/xg(A,R)+ / (g(A5(M),DTR)-g(DT(A5(n)),R)) 

where # denotes a contraction operation between tensors. The last term represents 
the contribution of the initial data on S. 



^•^These contractions result in a special structure of the quadratic term, crucial to the analysis 
in [Kl-Ro5] where we investigate a brealcdown criterion in General Relativity. The structure of 
this term is somewhat analogous to the corresponding term in the Yang-Mills theory, see previous 
section. 
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Representation (63) opens the possibility of proving a pointwise bound on the cur- 
vature tensor in terms of initial data on S and, as in the Yang-Mills case, the flux 
of curvature along the null boundary J\f~ [p) of the set VL. However, as opposed 
to the Yang-Mills equations on Minkowski background, where the curvature flux is 
bounded by the L^-norm of the curvature of initial data, no such a priori bounds 
are available for the Einstein vacuum equations. This suggests the use of the 
based curvature norms to deduce a breakdown criteria in General Relativity, i.e. 
to show that the space-time can be continued as long as such norms remain finite. 
For the Yang-Mills problem in Minkowski space the underlying reason for having 
an a- priori bounds on the flux of curvature is due to the presence of the Killing 
vectorfield dt = ■§^- In the case of the Yang-Mills equations on a smooth curved 
background, such as in [CS], the result remains true even though dt is no longer 
Killing; it suflices that its deformation tensor is bounded. We call such a vectorfleld 
approximately Killing. 

These considerations suggest the following question. Assume that the space-time 
(M, g) possesses an appivximately Killing, unit, vectorfield T, orthogonal to a 
space-like Cauchy hypersurface S, with deformation tensor 7r(X, Y) = g(DxT, Y). 
We also assume that the slices Sj obtained by following integral curves of T from 
So have constant mean curvature. Can the space-time be extended as long as tt 
remains finite in the uniform norm? 

The finitcnoss of the deformation tensor tt allows one to control, via energy estimates 
based on the Bel-Robinson tensor, both the norms of the curvature R along 
and the flux of curvature along the null boundaries M~{p). The key step in the 
remaining analysis is to derive a pointwise curvature based on the representation 
formula (63). In [Kl-Ro5] we give an affirmative answer to the question raised above 
by showing that the size of the region of validity of the formula (63) depends only 
on the assumed L°°- bounds on tt and reasonable assumptions on the initial data 
on E. Such estimates follow from our work in [Kl-Ro4]. More precisely we show 
that for a space-time metric g in the form 

g = —n^dt^ -\- gijdx''dx\ 

where n is the lapse function of the t foliation and the vectorfield T is orthogonal 
to St, the following result holds true. 

Theorem 5.2. Assume that (M,g) is a globally hyperbolic Einstein vacuum space- 
time with So a Cauchy hypersurface. Let the lapse function n and the deformation 
tensor irofT satisfy 

No^<n<No, IIttIIloc </Co 

Assume also that M contains a future, compact set I? C M such that for any point 
q the radius of injectivity of N~{q) is at least 5q> Q. 

Let S be one of the slices of the t foliation. Then assumptions Al, A2 of this paper 
are satisfied for all points p at distance < (5* from S, where (5* depends only on the 
Cauchy data on Sq, A^O; (ind 6o. In particular, the representation formula (63) 
holds for all such points. 
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6. Open questions 

All the results of this paper have been derived under assumption A2 which re- 
quires that, for any point p, the boundary A/"" (p) of the causal past J~{p) remains 
smooth at least until it reaches the space-like hypersurface E. It is only under this 
assumption that we can guarantee that the Kirchoff-Sobolev parametrix of Theo- 
rem 3.11 gives a faithful representation of a solution of the wave equation. We have 
already discussed the two obstructions to smoothness of M^{p): conjugate points 
of the congruence of past directed null geodesies from p and intersection of two 
distinct past directed null geodesies from p. The second obstruction can be easily 
demonstrated on a space-time M = M x or M = M x !!„ based on a flat torus 
or a flat cylinder Ha of width a. In those cases, however, an (exact) parametrix can 
be easily constructed by lifting the problem to the covering space R x IR^. On the 
other hand, the examples above are very special. Conjugate points can not be re- 
moved by a simple^** lifting and the quantity trx, which features prominently in our 
representation formulas, diverges to — oo at a conjugate point. However, the same 
focusing phenomenon shrinks the volume of a conjugate point region thus leaving 
open a possibility that, perhaps, with some additional assumptions on the structure 
and strength of conjugate points, the integral quantities appearing in our represen- 
tation formulas remain finite. It may thus be that the Kirchoff-Sobolev parametrix 
remains valid even beyond the region of formation of conjugate points. A good 
place to start investigating this issue would be product manifolds M = R x M 
with metrics of the form g = —dt^ + Qijdx^dxK A particularly interesting class to 
consider are product manifolds with M collapsing in the sense of Cheeger-Gromov. 
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